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Abstract 

Analytic solutions of the magnetohydrodynamic equilibrium equations for a cylindri- 
cally symmetric magnetically confined plasma with reversed magnetic shear, s < 0, and 
sheared flow are constructed by prescribing the safety factor-, poloidal velocity- and ax- 
ial velocity- profiles consistently with experimental ones. On the basis of the solutions 
obtained in most of the cases considered it turns out that an increase of \s\ and of the 
velocity components result in larger absolute values for the radial electric field, E r , its 
shear, \dE r /dr\ = \E' r \, and the E x B velocity shear, c^exb = \d/dr(E x B/5 2 )|, which 
may play a role in the formation of Internal Transport Barriers (ITBs) in tokamaks. In 
particular for a constant axial magnetic field, o;e x b at the point where E' r = is propor- 
tional to 1 — s. Also, \E' r \ and c^exb increase as the velocity shear takes larger values. The 
results clearly indicate that s < and sheared flow act synergetically in the formation of 
ITBs with the impact of the flow, in particular the poloidal one, being stronger than that 
of s < 0. 
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1. Introduction 



Understanding Internal Transport Barriers (ITBs) in plasmas is very important for 
the advanced tokamak scenarios The ITBs usually are associated with reversed 

magnetic shear profiles arid their main characteristics are steep pressure profiles in 

the barrier region jB] and radial electric fields associated with sheared flows [7j, jH]. The 
mechanism responsible for the formation of ITBs is far for completely understood. It is 
believed that the flow, the radial electric field, its shear and the E x B velocity shear, 



play a role in the barrier formation by mode decorrelation thus resulting in a reduction 
of the outward particle and energy transport [3], [TU] . 

The experimental evidence up to date has not made clear whether the reversed mag- 
netic shear, s < 0, or the sheared flow (toroidal or poloidal) are more important for the 
ITBs formation. In some experiments the safety factor profile is considered as the crucial 
quantity (e.g. JI]) while according to others the necessity of reversed magnetic shear 
is questionable (e.g. P). On the other hand, the flow-either toroidal [T2| or poloidal 
may be important in the formation of ITBs. Also, it has been argued that the 
toroidal velocity may be more important than the poloidal one (see for example Ref. [TH]). 
It should be noted, however, that only few direct measurements of the poloidal velocity 
have been performed; this velocity is usually calculated by means of neoclassical theory 



The aim of the present work is to contribute to the answer of the above mentioned open 
questions by studying magnetohydrodynamic (MHD) cylindrical equilibria with reversed 
magnetic shear and sheared flow. The study can be viewed as an extension of a previous 
one on tokamak equilibria with incompressible sheared flows and monotonically increasing 
q-profiles in connection with certain characteristics of the L-H transition [T^]. The work is 
conducted through the following steps: The profiles of certain free quantities, including the 
safety factor and the velocity components are first prescribed and then exact equilibrium 
solutions are constructed self consistently. This is the subject of Sec. 2. In Sec. 3 on the 
basis of the solutions obtained the equilibrium properties are examined and the impact 
of s < and the flow on E r , E' r and we x b is evaluated. The conclusions are summarized 
in Sec. 4. 

2. Cylindrical equilibria with reversed magnetic shear 

The equilibrium of a cylindrical plasma with flow satisfies (in convenient units) the 
relation 



stemming from the radial component of the force-balance equation p(v- V)v = j x B — VP 
with the aid of Ampere's law. Here, P is the plasma pressure; Bg and B z are the poloidal 
and axial components of the magnetic field, respectively; M| = (vgg)/Bg is the square of 
the Mach number defined as the ratio of the poloidal velocity to the poloidal-magnetic- 
field Alfven velocity. Because of the symmetry any equilibrium quantity depends only on 
the radial distance r and the axial velocity v z as well as the velocity shear do not appear 
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in (J2J); also, the flow is incompressible. In addition to v z four out of the five quantities in 
(J2J) can be prescribed. 

On account of typical experimental ITB profiles we prescribed the quantities q, B z , 
Vg, v z and g as follows: 

strongly reversed shear profile (SRS) (Fig. |TJ) 

q{p) =q c (l — p + — p (3) 

V "c r m i n qc r m m J 

or alternatively weakly reversed shear profile (WRS) 

q{p) =qA\ p + — p (4) 

V (fc r min qc ?"mm / 

where p = r/R Q with r defining the plasma surface, q c = q(r = 0), r min is the position 
of minimum q, and Aq = q c — g min . The SRS-profile (jSJ) does exhibit a maximum at 
the plasma center r = in addition to the minimum one at r = r m ; n and has stronger 
magnetic shear in the central region just inside the q m { n position than that of the WRS 
one. It should be clarified, however, that the WRS profile (@J), which does not have an 
extremum on the magnetic axis r = 0, has been chosen in order to simplify the calculations 
though the physical situation may not be well represented in the immediate vicinity of 
the magnetic axis; 

B z = B z0 [l + 5(l-p 2 )} 1/2 (5) 

where B z0 is the vacuum magnetic field and the parameter 5 is related to the magnetic 
properties of the plasma, i.e. for 5 < the plasma is diamagnetic; 
Gaussian-like poloidal velocity profile 



v e = 4v 90 p(l - p) exp 



(p Pmin) 
h 



(6) 



where the parameter h determines its broadness and vgo is the maximum of vg; either 
peaked axial velocity profile 

^ = M1-P 3 ) 3 (7) 
or Gaussian-like v z profile similar to that of (jHJ); and the density profile 

£? = £?o(1-P 3 ) 3 - (8) 

The following quantities can then be calculated: the poloidal magnetic field Bq = 
epB z /q where e = t /Rq is the inverse aspect ratio with 2ttR associated with the length 
of the plasma column; the magnetic shear s = (r/q)(dq/dr)\ the current density via 
Ampere' s law; the electric field via Ohm' s law; its shear E' r and cuexb 

by (HJ). Also, 

integration of (0) so that P(r = r ) = yields the pressure. The calculations have been 
performed analytically by developing a programme for symbolic computations ^U] m 
connection with Ref. ^7j. This also allowed us to examine conveniently purely poloidal 
flows, purely axial flows, 2-pinch configurations or ^-pinch configurations as particular 
cases. The analytic expressions which can be derived readily by the programme are 
generally lengthy and will not be given explicitly here. Some concise and instructive 
expressions will only be presented in the next section along with typical profiles for the 
calculated quantities supporting the results obtained. 
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3. Results 



We have set the following values for some of the parameters: B z0 = 1 Tesla, Qq = 8.35 x 
10~ 8 kgr/m 3 corresponding to n = 5 x 10 19 particles/m 3 , p min = 0.5, e = r /_R ~ 1/3? 
5 = —0.0975, g m i n = 2, max^ = 1 x 10 4 m/sec and maxv z = 1 x 10 5 m/sec; Consequently, 
it is guaranteed that M| « M|, where M| = (v^gj/B^, a scaling typical in tokamaks 
because B z « 10 and f 2 w 10 i^. It is noted here that since in tokamaks Mg < 0.1 the 
flow term in (J2J) is perturbative around the "static" equilibrium Mg = 0. Also, the choice 
<?min = 2 was made because according to experimental evidence for g min < 2 strong MHD 
activity destroys confinement possibly due to a double tearing mode [TH] . A similar result 
was found numerically for one-dimensional cylindrical equilibria with hollow currents in 
Ref . ^H] ■ The impact of the magnetic shear and flow on the equilibrium, in particular on 
the quantities E r , E' r and uj ExB , was examined by varying the parameters q c , Aq, h, v z0 , 
and v eo [Eqs. ©, Q, © and ©]. 

For reversed magnetic shear profiles we came to the following conclusions: 



1. Pressure 



Substitution of B e and its derivative in terms of q and s in (j2J yields 

p '- B ^[ 1+ ( £ i) 2 ] +ropK2+(s - 2)] (^) 2 ' (9) 

For s < 0, increase of \s\ makes the pressure profile steeper (see also Fig. 2). Equa- 
tion (jnj) also implies that the pressure profile becomes steeper when the plasma is more 
diamagnetic, i.e. when B' z related to the parameter 5 in (JHJ) takes larger values. 

2. Current density 

• The axial current density profile becomes hollow and, irrespective of the reversal of 
the magnetic shear, there is a critical distance p cr outside the g m i n position at which 
J z becomes negative (Fig. 3). In particular, for B z = B z $ = const, one obtains 

J, = ~{tB 9 ) = ^{2-s) (10) 
r ar R^q 

Consequently, for s > 2, J z reverses. The radial distances at which J z = for the 
SRS [Eq. Q] and the WRS [Eq. (JJJ)] q-profiles, respectively, are 



1/3 

SRS 
rcr 



and 

WRS = . Vc_ 
Her Pmin ^ • 

Therefore, the position of p cr is shifted towards the center as s takes lower nega- 
tive values. It is noted here that equilibrium toroidal current density reversal for 
monotonically increasing g-profiles was reported in Ref. [20] (Fig. 3 therein). 



4 



Very large values of Aq on the order of 10 2 result in the formation of j z profiles with 
"holes" in the central region- j z m 0- inside the p m \ n position as demonstrated in 
Fig. 4, a result consistent with experimental evidence 



The total axial current I z = 27rr Bg(r ) for SRS profiles is smaller than that for 
WRS profiles. 



3. E r and EL 



r 



Typical E r profiles exhibit an extremum in the region around q m in arid vanish at 
p = and p = 1 in agreement with experimental ones [7j, j2H|. Profiles with more 
than one extrema are also possible in the case of peaked v z profiles, localized vg ones 
and v z vq > as demonstrated in Fig. ED Experimental profiles of this kind were 
reported in Ref. (Fig. (9) therein). 

The main contribution to E r comes from the velocity, to which is proportional, and 
particularly from the poloidal one (Fig. |UJ). 

E r is sensitive to the relative orientation of v z , Vg and B z \ in particular, for v z vg < 
\E r \ is larger than that for v z vg > 0. (Fig. EJ). Similar results hold for E' r (Fig. EJ) 
and cjexb • 

For extended velocity profiles with v z ^ 0, an increase of \s\ results in an increase 
of \E r \ (Fig. ED, \E' r \ and wexB' If v z — 0, however, |s| has no impact on \E r \ and 
| |, as can be seen by inspection of E = v x B, and very weak impact on c^exb- 
This result indicates that the presence of v z "activates" the impact of s on E r , E' r 
and w E xB. 

An increase of the velocity shear nearly does not affect or even decreases the maxi- 
mum \E r \ (Fig. ITUj) but increases \E' r \ (Fig. [TTj) . 



4- CJexB 



A typical profile of c^exb has two large local maxima at the positions where the 
edges of the barrier are expected to be located in addition to other two smaller local 
ones (Fig. [T2j) . In most of the cases considered the maximum in the s < region 
is slightly larger than that in the s > region, (see Fig. IT2]) . In particular, for 
B z = const, at the point where E' r = one obtains: 



w ExB 



S)(^-Vg) 



R q 



:n) 



Eq. implies the following: 



1- wexb depends on the relative sign of v z , vg and B z , a result which we confirmed 
by wexb profiles obtained via the symbolic computation programme. 

2. The factor (1 — s) indicates that wexb for nearly shearless stellarator equilibria 
may be lower than that for tokamak equilibria with s < 0. 
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3. Despite the scaling v z « 10vg, for tokamak pertinent parametric values the 
contributions of v z -in connection with the term epv z /q- and vq to uexB are of 
the same order of magnitude, a result indicating the importance of the poloidal 
velocity. 

• For extended velocities (large values of the parameter h or/and peaked v z - profile) 
a percentage increase of \s\ in the barrier region results: 

1. in approximately the same percentage increase of we x b if the velocity is purely 
axial (Fig. E3J). 

2. nearly does not affect the value of wexb if v$ ^ 0. 

• An increase of the flow shear (variation of the parameter h from 0.1 to 0.001) causes 
a mean percentage increase of c^exb as large as 0.7 of that of the flow shear. (Fig. 
d. 

• The impact of a variation of the vg-shear on WexB is stronger than of the same 
variation of the t> z -shear. 

• The maximum increase of we x b is caused by vg in the case of non- vanishing peaked 
v z profiles. 

• Inspection of ve x b =ExB/5 2 and (0) implies that uexb for a z-pinch is equal 
to that for an equilibrium with purely axial flow. The same equality is valid for a 
#-pinch in comparison with an equilibrium with purely poloidal flow. In addition, 
it holds that 

^ExB-z-pinch ~ lO^ExB-e-pinclr (^) 
4. Conclusions 

The self consistent study of cylindrical equilibria with reversed magnetic shear and 
sheared flow presented in the previous sections led to the following conclusions: 

1. For reversed magnetic shear profiles (s < 0): 

• The larger values of \s\ the steeper the pressure profile. 

• The axial current density profile become hollow. 

• Strong reversed shear profiles formed by appropriately large values of Aq are 
associated with "hole" axial current density profiles. 

These results are consistent with experimental ones. 

2. Irrespective of the sign of s the axial current density can reverse in the outer plasma 
region, the reversal point being shifted towards the plasma core as s takes lower 
negative values. 

3. An increase of either \s\ or the velocity results generally in an increase of \E r \, \E' r \ 
and co>exb- 

4. An increase of the velocity shear results in an increase of \E' r \ and c^exb- 
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5. For a given value of \s\, c^exb takes slightly larger values in the s < region than 
in the s > region. 

6. E r , E' r and cjexb are sensitive to the relative orientation of vg, v z and B z . In 
particular, they take larger values for v z vg < rather than for v z vg > 0. 

7. The presence of v z activates s < 0, in the sense that for v z = 0, E r and E' r are 
s-independent. Also, for v z = 0, s < has very weak impact on u>exb- 

8. The impact of the poloidal flow and its shear on E r , E' r and cj ExB is stronger than 
that of the axial flow and the magnetic shear. 

Presuming that E r , E' r and cj ExB are of relevance to the ITBs formation, the above 
results clearly indicate that the reversed magnetic shear and the sheared flow have syner- 
getic effects on this formation with the flow, in particular the poloidal one, and its shear 
playing an important role. 
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Figure captions 



Fig. [TJ SRS and WRS safety factor profiles associated with Eqs. (J3J) and (jU), respectively. 
It is noted that the finite slope of the WRS curve at p = may not represent well the 
physical situation in the immediate vicinity of the magnetic axis. 

Fig. El WRS pressure profiles for Ag = 4 and Ag = 14. 

Fig. El Toroidal current density profiles for Aq = 4. It is noted that the finite slope of 
the WRS curve at p = may not represent well the physical situation in the immediate 
vicinity of the magnetic axis. 

Fig. HI Toroidal current density profile for WRS, q c = 102 and Aq = 100 that demon- 
strates the current "hole' in the core region. 

Fig. 03 Electric field profile for WRS with v z peaked and v g localized having three local 
extrema. 

Fig. IH1 Two -EV-profiles the one with v z — and the other with v g = for SRS and 
Gaussian-like velocity profiles. E r is normalized with respect to its value at p = 0.5 for 
v z = . 

Fig. [7| Two -E r -profiles for v z peaked and SRS, the one with v e • v z > and the other 
with vg ■ v z < 0. The profiles are normalized with respect to the first case at p = 0.5. 

Fig. |S1 Two profiles of E' r with v z peaked and SRS, the one with vg -v z > and the other 
with v g ■ v z < 0. The profiles are normalized with respect to the first case at p = 0.3. 

Fig. El Profiles of E r with peaked axial and extended poloidal velocities for WRS and 
two different values of Ag. The profiles are normalized with respect to the case with 
Aq = 4 at p = 0.5. 

Fig. El Profiles of E r with v z — for SRS and either extended (h = 0.1) or localized 
(h = 0.001) poloidal velocity. The profiles are normalized with respect to the first case at 
p = 0.5. 

Fig. ^3 Two profiles of E' r for SRS with v z — the one for extended (h = 0.1) and 
the other for localized (h = 0.001) poloidal velocities. The profiles are normalized with 
respect to the second case at p = 0.55. 

Fig. H21 Typical '-'-'ExB'Profile for WRS, peaked axial and localized poloidal velocities. 

Fig. [131 Profiles of c^exb for WRS, peaked axial velocity, and either Ag = 4 or Ag = 14. 
The profiles are normalized with respect to the first case at p = 0.45. 

Fig. EH c^ExB-profile for SRS, Gaussian-like axial and poloidal velocity components 
both either extended (h = 0.1) or localized {h = 0.001). The profiles are normalized with 
respect to the first case at p = 0.3. 
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Figure 10: 
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